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Abstract: A 891 g BGO (Bi 4 Ge 3 0i2) scintillating bolometer has been tested at 10 mK in the 
underground Laboratori Nazionali del Gran Sasso (Italy). The discrimination capability, the radio- 
purity of the compound and the main features of the crystal have been studied in order to demon- 
strate the excellent performances obtained by operating a scintillating bolometer in the field of y 
and a spectroscopy. The sensitivity of this detector in the study of extremely low surface contami- 
nations has been investigated. 
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1. Introduction 

The neutrinoless double beta decay (OvDBD) has a widely recognized importance in modern 
physics, as its direct observation would demonstrate the Majorana nature of neutrino and set a 
value on the absolute mass of this particle. 

Cuoricino, designed for the purpose of studying the OvDBD of 130 Te, was operated from 2003 
to 2008 in the Laboratori Nazionali del Gran Sasso (Italy) and is still one of the most sensitive 
experiments ever realized in this field UU] . Cuoricino was based on the bolometric technique and 
consisted of 62 Te02 crystals (corresponding to 1 1.3 kg of 130 Te) operated at a constant temperature 
of ~10 mK. Thanks to the excellent performances provided by this detector in terms of efficiency 
and energy resolution, Cuoricino could set an upper limit on the neutrino Majorana mass of 300 - 
710 meV, depending on the adopted nuclear matrix element evaluation. 

The important physics results provided by this experiment lead to its upgrade, CUORE, that 
will improve the discovery potential on the effective Majorana neutrino mass to 41 - 95 meV [Q, ED- 
The achievement of this goal is subject to the suppression of the background counting rate in the 
energy region of the 130 Te OvDBD (~2527 keV) down to 10~ 2 counts/keV/kg/y, which means 
about an order of magnitude less with respect to Cuoricino. 

Such a low background can not be reached only by locating the experiment in deep under- 
ground laboratories and by equipping it with proper shields. The Cuoricino experience, indeed, 
showed that the most important sources of background are the radioactive contaminations of the 
materials located close to the detector. For this reason, the materials of CUORE must have a bulk 
contamination in 238 U e 232 Th (or their daughters) lower than a few pg/g, as well as surface con- 
taminations lower than a few nBq/cm 2 . Unfortunately, due to the low level of the contaminations, 
these sources of background are very difficult to identify and reduce. 

Bulk contaminations of the materials can be studied with a good sensitivity through conven- 
tional techniques, i.e. high purity germanium (HpGe), inductively coupled plasma mass spectrom- 
etry (ICP-MS) and neutron activation. 
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The situation is much more complicated for the material surface contaminations. Since "stan- 
dard" techniques (i.e. Si Surface Barrier Detectors) can not achieve a sufficient sensitivity, the only 
way to access these extremely low contaminations is the bolometric technique. Bolometers, in- 
deed, have no dead layer and can reach energy resolutions orders of magnitude better than standard 
Silicon Surface Barrier Detectors. 

Unfortunately, pure bolometers can not identify the nature of the particles that interact in the 
detector and so the characterization of the background sources remains a very difficult task. 

In this paper we demonstrate how this problem can be overcome by means of a BGO scintil- 
lating bolometer. 

2. Scintillating Bolometers 

Bolometers are calorimeters in which the energy released in the absorber by an interacting particle 
is converted into phonons and measured via a temperature variation 0]. An ideal instantaneous 
deposition of energy E in the absorber gives rise to a temperature variation AT oc E/C, where C is 
the thermal capacitance of the bolometer itself. A small capacitance is therefore needed to obtain 
sizable AT. For this reason, bolometers are usually made of dielectric and diamagnetic crystals and 
are operated at cryogenic temperature (»10 mK). 

The temperature increase is measured through dedicated sensors, whose resistance shows a 
steep dependence on the temperature. For Neutron Transmutation Doped (NTD) Ge thermistors 
05D, the resistance variation is than converted into a readable voltage signal by biasing the thermis- 
tors with a constant current. 

If the bolometer is a scintillating crystal, the heat signal can be combined with the light sig- 
nal. The composite device shows many advantages: the bolometric technique provides an excellent 
energy resolution and a high efficiency, while the simultaneous read-out of the light allows to iden- 
tify the interacting particles (n, a, jS/y/jU) thanks to their different Light Yield (LY). Indeed, the 
bolometers are very slow detectors (pulse time development of the order of hundreds of millisec- 
onds) and do not allow to recognize the interacting particle by the shape of the thermal signal except 
in some cases in which the energy release in the crystals produces states with very long decay times 

01. 

Examples of other interesting features of scintillating bolometers can be found in [J7fl. 
Our purpose is to demonstrate that an accurate characterization of the surface contaminations of a 
material can be obtained in a simple way by facing the material itself to a scintillating bolometer. 
The precise measurement of the particle energy, as well as the tagging of the particle nature, allows 
for the identification of the single contaminants. At the same time, the large active surface and the 
absence of dead layers provide a very high sensitivity. 

In this paper we present the results of a measurement performed with a 5x5x5 cm 3 BGO 
scintillating crystal (891 g). The choice of this crystal is due to many interesting features: first of 
all, the radio-purity of the compound and the large light yield. Then, the high atomic number of 
BGO and the large dimensions of our sample, which make this detector very appealing also for y 
spectroscopy. 
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3. Detector 



The detector we present is a BGO crystal whose surfaces are polished at optical grade. The crystal 
is held by means of L-shaped Teflon pieces (PTFE) fixed to copper columns, and it is surrounded 
by a reflecting foil (3M VM2002) in order to increase the light collection efficiency (see figure [I]). 




Figure 1. The 5x5x5 cm 3 BGO crystal surrounded by the reflecting foil. 

The emitted light is measured by a second bolometer. This Light Detector (LD) is realized 
with a 36 mm diameter and 0.3 mm thick pure germanium disk covered with a 600 A thick layer of 
Si02 in order to ensure a good light absorption. 

The composite detector is mounted in the Oxford 200 3 He/ 4 He dilution refrigerator located in 
the Laboratori Nazionali del Gran Sasso («3650 m w.e.). The details of the cryogenic facility can 
be found in OSfl. 

The signals produced by the thermistors are read via a cold pre-amplifier stage located inside 
the cryostat. A second stage of amplification, as well as the front end, are located on the top of 
the cryostat at room temperature. The signals are filtered by means of an anti-aliasing Bessel filter 
with cut frequencies of 63 Hz for the BGO signals and 120 Hz for the Light signals [@J. The entire 
waveform is digitized and recorded with a sampling frequency of 2 kHz. 

The trigger of the BGO and of the LD is software generated. The LD is acquired in coincidence 
with the heat channel in order to reject most of the signals that are not due to a BGO scintillation 
event. 

The off-line analysis allows to determine the pulse amplitude, as well as many pulse parame- 
ters, by means of the optimum filter technique HTQl [TTTl, 
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4. Results 



The BGO crystal cooled very slowly. A similar behaviour was noted also in a previous measure- 
ment with 2x2x2 cm 3 BGO crystals. This behaviour can be ascribed to some excess heat capacity 
that decouples at low temperatures. In spite of the cooling problem, it was possible to perform an 
accurate measurement with this crystal. 

The light vs heat scatter plot obtained in w455 hours of measurement is reported in figure 0. 
The excellent discrimination between a and /} /y events can be observed. 

The heat spectrum (reported in figure ||) is calibrated on the 208 T1 line and on the y peaks of 207 Bi 




Energy (keV) 

Figure 2. Light vs heat scatter plot of the 5x5x5 cm 3 BGO crystal for a w455 hours background measure- 
ment. In blue /3/y events, in orange a events. 

which is produced by proton interactions on 206 Pb HT2H . It is probably mainly originated by the 
interaction of cosmic ray protons since the BGO crystal was produced in China (SICCAS) and 
delivered by plane. 

The light spectrum (figure is calibrated by means of a 55 Fe source faced to the LD. The 
calibration is performed with a simple linear function, by attributing the nominal energy of the 
55 Fe X-rays (5.89 keV and 6.49 keV) to the identified peaks. 

The radio-purity of the crystal has been studied by investigating its internal contaminations. 

In the a region several peaks can be identified. These are due to internal contaminations of 

209 Bi> 210m Bi and 210p Q _ 

2i0mgj j s p roc j uce d v j a the 209 Bi(n,y) 210m Bi reaction and accumulates in the crystal due to its 
long half-life (T 1/2 =3 x 10 6 y) Q. 210m Bi alpha-decays on 206 T1 with a Q-value of «5.3 MeV. The 
decay can be on the fundamental level, giving rise to the lowest spot at 5.3 MeV (see figure @), or 
on the excited state of 206 T1 with the prompt emission of y rays of ~300 and ~640 keV. Since the 
probability that these y rays release all their energy in the BGO crystal itself is high (>54%), these 
mixed a+y events give rise to other two spots that lie between the a and the /3/yband. 

209 Bi (Q=3137 keV) was supposed to be the heaviest stable isotope, until P. de Marcillac et. 
al detected its decay HI4H . Thanks to the simultaneous observation of the ground and excited state 
decay of this isotope, our measurement provides the ultimate test on the rare decay of 209 Bi IH3I1 . 
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Figure 3. Heat Spectrum of a 5 x 5 x 5 cm 3 BGO crystal for a «455 hours background measurement. Blue: 
j3/y events. Orange: a events. 
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Figure 4. Light Spectrum of a 5 x 5 x 5 cm 3 BGO crystal for a «455 hours background measurement. Blue: 
total spectrum. Orange: a events. 

Finally, the high number of counts due to the decay of 210 Po (Q=5407 keV) is probably due 
to the activation of 209 Bi by the 209 Bi(n,y) 210 Bi reaction. 210 Bi subsequently gives origin to 210 Po 
through j6-decay (Ti/ 2 ~ 5 days). 

The full understanding of these contaminants allows for a detailed analysis of the resolution 
and of the light emission of the crystal. 

The FWHM energy resolution of the heat channel is (37.5 ± 0.5) keV, with no dependence on 
the energy of the peak. This value is obtained using the energy spectrum of the entire measurement. 
However, better results can be obtained fitting only the spectrum of the last days of measurement, 
as the cooling of the crystal determined an increase in the signal to noise ratio and, therefore, an 
improvement in the energy resolution. In a couple of weeks of measurement, indeed, we observed 
an in improvement of the resolution of 51%. 

The performance of the light detector is excellent. The energy FWHM resolution of the light 
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Band 


Energy 


Background 




[MeV] 


[c/keV/kg/y] [nBq/cm 2 /MeV] 


P 


3.0 - 4.0 


0.01 - 0.07 — 


a 


2.0-3.0 


0.03-0.14 6.32-30.45 


a 


3.5-4.0 


0.04 - 0.22 8.50 - 48.84 



Table 1. Background obtained with the 5x5x5 cm 3 BGO crystal in three different regions. The values are 
evaluated with Feldman-Cousins at 68.27% C.L. The results were obtained by assuming that all the a events 
were due to surface contaminations. 

signals ranges from 0.78 keV to 1.5 keV (the first value calculated on the 570 keV peak of 207 Bi, 
the second one on the 2615 keV line of 208 T1). 

The analysis of the light emission of the BGO crystal demonstrates the discrimination capa- 
bility of this detector. 

The LY of the 7 peaks is evaluated on 40 K, 208 T1 and each line of 207 Bi and no dependence on 
the energy is observed. The average value is LY(j8/y) = (16.61 ± 0.02) keV/MeV. 

On the contrary, the LY of a particles increases with the energy: we measured LY = (2.482 
± 0.002) keV/MeV at 3137 keV ( 2()9 Bi line) and LY = (3.011 ± 0.003) keV/MeV at 5400 keV 
( 210 Po line), corresponding to a quenching factor for a particles of 0.1494 ± 0.0002 and 0.1813 
± 0.0002, respectively. This behavior, which is a consequence of the energy dependence of the 
stopping power (Birks law HTBH ). was already observed in other crystals HT7H . 

The low values of the quenching factor for a's, together with the narrow bands allow to dis- 
criminate very well a and /$/y events. 

In addition, except for the 207 Bi, 209 Bi and 210m Bi isotopes, no evidence of internal radioactive 
contamination of the crystal is observed. In particular natural radioactive contaminations have an 
activity lower than 0.14 pg/g in 238 U and 0.41 pg/g in 232 Th at 90% C.L. Most of the detected 
contaminations can be easily reduced by avoiding the shipping of the crystals by plane; in this way 
the cosmic ray activation can be strongly suppressed. 

The measurement presented in this paper allowed to perform a preliminary study of the sur- 
face contaminations. The analysis was focused on the a band in the 2-3 MeV region and in the 
3.5-4 MeV region (the division is due to the presence of the 209 Bi peak). The study of surface con- 
taminations can be of particular interest for OvDBD experiments, whose goal is the achievement 
of a low background counting rate in this energy region. The results are reported in Table [I]. The 
evaluation of the surface activity is obtained assuming conservatively that all the alpha events in 
the selected energy regions are due to surface contaminations of materials facing the detector. The 
background obtained in the jS/y region with E>3 MeV is also reported. The measured value also 
includes events due to pile-up of /$/y due mainly to the high internal contamination of 207 Bi. As 
mentioned above this rate can be reduced simply by avoiding the shipping of the crystals by plane. 

The j3/y region with E>3 MeV is very interesting especially for experiments aiming to study 
the OvDBD of isotopes with large Q-value. Indeed, above this line the environmental background 
due to natural y's decrease abruptly. 

The very low counting rate obtained in the energy region of interest is rather surprising if one 
considers that we did not apply any cleaning procedure on the crystal and on the reflecting foil that 
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was surrounding it. 

Since this technique will be applied to the study of surface contaminations, we performed 
Monte Carlo simulations (GEANT-4) to understand the efficiency of this detector and then its 
sensitivity. The simulations were performed with an array of four 5x5x5 cm 3 BGO crystals faced 
to a slab of material to be analyzed (figure [5]). The goal of these simulations is to study the geometric 
efficiency (e), i.e. the ratio between the number of alpha particles that hit the detectors and the 
number of generated particles on the surface facing the detectors. For this purpose, the simulations 
were performed generating a decays on the surface of the material (no depth). 

In addition, different simulations were performed with a particles with energy between 3 and 
6 MeV to study the dependence of the geometric efficiency on the particle energy. The simulations 
confirmed our expectations, showing that the energy dependence is less than 1% (limit given by 
the simulation statistics). 




Figure 5. Simulated array of four 5x5x5 cm 3 BGO crystals faced to a slab of material positioned under the 
crystals. 

The geometric efficiency for the detection of surface contaminations of a circular/square slab 
with a diameter/side of 12.4/1 1 cm faced to the BGO crystal array at a distance of 5 mm was evalu- 
ated to be e ~35%. For a measurements of 1 month with a sample of 121 cm 2 faced to the crystals 
this allows to achieve sensitivity of 21 nBq/cm 2 . The obtained sensitivity decreases with increasing 
depth of the contamination because of the worsening of the detection efficiency. 

In order to compare the sensitivity that could be reached with BGO crystals with respect to a 
standard a detector (i.e. Si surface barrier detector) and pure thermal bolometers (Te02 bolome- 
ters), several simulations were performed. The geometry of the simulation of the Te02 crystals 
was the same used for the simulation of the BGO. For what concerns surface barrier detector, we 
simulated a 4 cm diameter Si disk flTPH faced to the same slab of material to be analyzed. 

The results are reported in figure || for surface contaminations in 232 Th and 238 U with an 
exponential density profile and a depth A of 5 micron (i.e. the one used for Cuoricino background 
reconstruction [EOT ). 
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Figure 6. Simulations of energy spectra that can be obtained with a surface contamination in 232 Th (on the 
top) and 238 U (in the bottom) with an exponential density profile and a depth A of 5 micron by a Si surface 
barrier detector (green), an array of four TeC>2 crystals (blue) and an array of four BGO crystals. For the last, 
thanks to the possibility to identify the interacting particle, the beta/gamma region (orange) and the alpha 
region (red) are reported separately. 

As can be seen in the figures, the main advantage of the BGO detector is the possibility to per- 
form separately a and j3/y spectroscopy. This feature allows to study with a very high sensitivity 
low energy a's, that in the case of Te02 and Si detector would be overwhelmed by the y natural 
radioactivity. At the same time, with the BGO detector it is possible to access rare y's and /3 events 
in the energy region dominated by the a background (above 2.6 MeV). 

Another important advantage of the BGO detector is the very low intrinsic background, which 
is about two orders of magnitude lower with respect to the count rate of a standard Si detector (0.05 
counts/h/cm 2 between 3 and 8 MeV lUPIl ). This leads to a much better signal to noise ratio and, 
therefore, to a much larger sensitivity on surface contaminations. 

Finally it must be stressed that, unlike Si detectors, bolometers can easily achieve large ac- 
tive surfaces and do not have dead layers. Moreover, the excellent energy resolution provided by 
bolometers (~5 keV FWHM energy resolution with a Te02 bolometer flZTH with respect to 25- 
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30 keV FWHM resolution of a Si detector) makes these detectors particularly appealing for high 
sensitivity spectroscopy. 

5. Conclusion 

For the first time a large BGO crystal (891 g) was measured at cryogenic temperature, giving 
excellent results in terms of discrimination capability and radio-purity. 

We demonstrated that this crystal is a perfect candidate to perform y and a spectroscopy with 
a very high sensitivity. 

Looking at these results, we are planning a measurement with an array of 4 BGO crystals 
(5x5x5 cm 3 ). The array will allow for a coincidence analysis and, therefore, for a precise de- 
termination of the surface contamination of the crystals. In addition the reflecting foil will be 
removed. In this way it will be possible to measure the surface contaminations of the materials 
close to the detector with a sensitivity of « nBq/cm 2 in few weeks. 
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